matlab code bisection method